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Abstract 

The ground-state hyperfine spHtting values of high-Z boronhke ions are calculated. Calculation of the 
interelectronic-interaction contribution is based on a combination of the 1/Z perturbation theory and the 
large-scale configuration-interaction Dirac-Fock-Sturm method. The screened QED connections ai^e evalu- 
ated utilizing an effective screening potential approach. Total hyperfine splitting energies are presented for 
several B-like ions of particular interest: ^^Sc^^+, ^''Fe^^+, ^°''Pb^^+, and ^°^Bi^^+. For lead and bismuth 
the experimental values of the Is hyperfine splitting are employed to improve significantly the theoretical 
results by reducing the uncertainty due to the nuclear effects. 

PACS numbers: 31.30.Jv, 31.30.Gs, 12.20.Ds 
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I. INTRODUCTION 



First accurate calculations UJ, |2|] of the hyperfine splitting in highly charged ions were stimu- 
lated by astronomical search in hot astrophysical plasma ||3|]. Later, high-precision measurements 
of the ground-state hyperfine splitting in heavy H-like ions were performed [|4i|5i|6iLZil8|]. The main 
goal of such experiments was to probe quantum electrodynamics (QED) in the strong Coulomb 
field induced by a heavy nucleus. However, accurate theoretical calculations (see Ref. [9] and 
references therein) showed that the uncertainty of the theoretical results, which mainly originates 
from the nuclear magnetization distribution correction (the Bohr-Weisskopf effect), is comparable 
with the QED correction. For this reason, any identification of QED effects on the hyperfine split- 
ting in heavy H-like ions turned out to be unfeasible. It was shown, however, that this uncertainty 
can be significantly reduced in a specific difference of the hyperfine splitting values of H- and 
Li-like ions with the same nucleus lioll . High-precision measurements of the hyperfine splitting 
in heavy Li-like ions are presently in preparation llll[l . 

The motivation for accurate calculations of the hyperfine sphtting in B-like ions is twofold. 
From one side, high-precision prediction of the hyperfine splitting of B-like Fe may be important 
for astronomical search [|l2n . From the other side, the study of the hyperfine splitting in heavy 
B-like ions can be used to reduce the uncertainty associated with the Bohr-Weisskopf effect in 
some specific difference of the hyperfine splitting values for B- and Li-like ions or B- and H-like 
ions. The origin of this reduction is essentially the same as for the related (yf-factor values [1^ and 
can be easily seen from the approximate analytical expressions for the Bohr-Weisskopf correction 



given in Refs. 11141 . 



In this article we calculate the ground-state hyperfine splitting of B-like ions. The 
interelectronic-interaction correction of the first order in 1 /Z is evaluated within the rigorous QED 
approach. The higher-order terms are calculated within the large-scale configuration-interaction 
Dirac-Fock-Sturm method. QED corrections are calculated using an effective potential approach 
in order to account for the effect of screening. The experimental values of the Is hyperfine split- 
ting in H-like 207p|^8i+ 209gj82+ employed to evaluate the Bohr-Weisskopf correction for 
the corresponding B-like ions. Due to the correllation between the values of this correction for Is 
and 2pi/2 states mentioned above, the uncertainties of the theoretical values caused by the nuclear 
effects are strongly reduced. 

Relativistic units (h = c = 1) and the Heaviside charge unit (a = e^/(47r), e < 0) are used 
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throughout the paper. 



n. BASIC FORMULAS AND CALCULATIONS 



Within the point-dipole approximation, the hyperfine interaction is described by the Fermi- 
Breit operator 

I e I (cK ■ [/X X r]) 



Atx 



(1) 



where the vector a incorporates the Dirac matrices and /x is the nuclear magnetic moment operator. 

The ground-state hyperfine splitting of a B-like ion can be written in the following form 1 15, 1^ 

aiaZY m a 21 + 1 9 
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where m is the electron mass, nip is the proton mass, fi^ is, the nuclear magneton, and / is the 
nuclear spin. A^{aZ) is the one-electron relativistic factor 

6[2(l+7)- ^2(1 + 7) 
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7)^7(47^ 



(3) 



where 7 = a/1 — {aZy, 6 is the nuclear-charge-distribution correction, and e is the one-electron 
nuclear-magnetization-distribution correction (the Bohr-Weisskopf effect). The terms B^{aZ)/Z 
and C^{Z, aZ)/Z'^ determine the interelectronic-interaction correction to the first and higher or- 
ders in 1/Z, respectively. The Xrad term stands for the QED correction. 



Finite-nuclear-size correction 6 can be calculated both analytically lll4l ll7h and numerically. In 



the present work, it is obtained numerically by solving the Dirac equation with the Fermi model for 
the nuclear charge distribution. The Bohr-Weisskopf correction e is calculated within the single- 



particle model as described in Refs. , uM. Apart from this direct evaluation, at the end of 
this section we also derive the e values for 207p|^77+ ^j^^ 209gj78+ y^^[^g (-j^g experimental results 
for the Is hyperfine splitting. The nuclear root-mean-square radii are taken from Ref. and the 
nuclear magnetic moments are taken from Ref. ll20n . 

The QED correction to the hyperfine splitting of the first order in a consists of two parts, self- 
energy correction and vacuum-polarization correction. The self-energy correction (Fig. [T]) is the 
sum of irreducible, reducible, and vertex parts. 
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Figure 1 : Feynman diagrams representing the self-energy correction to the hyperfine splitting. The double 
line indicates the bound-electron propagator and the dashed line terminated with the triangle denotes the 
hyperfine interaction. 
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Figure 2: Feynman diagrams representing the vacuum-polarization correction to the hyperfine splitting. 
Notations are the same as in Fig.[T] 

where C^Mrjm '^he Clebsch-Gordan coefficients and ximi is the nuclear wave function. The 
irreducible part is given by the expression 



{a'\iT.{ea)--1%m)\n){n\H^,\a) 



(5) 
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Here \a) and \a') are the states of the valence electron with angular momentum projections m and 
m', respectively, is the energy of the state \n) in the binding potential under consideration, iJ^ 
is the magnetic-dipole hyperfine-interaction operator ([T]), and 6m is the mass counterterm. T.(e) 
denotes the unrenormalized self-energy operator defined as 



(a|S(e)|6) 
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where I{uj, xi, X2) = e'^a^a'^D^^{LU, X1 — X2) with the Dirac matrices = (1, a), and the photon 
propagator D^jj,. To separate the ultraviolet divergencies, the expression ([5]) is decomposed into 
zero-, one-, and many-potential terms. The zero-potential M^^J and one-potential M^^J terms are 



calculated in momentum space using the formulas from Ref. n2u\ . The residual part of Mi^i, the 



so-called many-potential term M^^^'^\ is calculated in coordinate space. The expressions for the 



reducible and the vertex parts read 

Mred = {a'\H^\a){a\-^J:{e) 
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Both reducible and vertex parts are ultraviolet-divergent, whereas the sum M„ = Mred + M^er is 



finite. Following Refs. 11221 . 



23L 



24, 



2511 . we separate out zero-potential term M^r^ and evaluate it 



in momentum space. The remaining many-potential term mH^^ is calculated in coordinate space 
as a point-by-point difference between the contributions with bound and free propagators in the 
self-energy loop. 

The angular integration and the summation over intermediate angular projections in the many- 
potential terms M^^^^^ and mH^^ are carried out in a standard manner. The many-potential terms 
involve infinite summation over relativistic angular quantum number k = ±(j + 1/2). The sum- 
mation is terminated at a maximum value \ k\ = 10, while the residual part of the sum is evaluated 
by the least-square inverse-polynomial fitting. For any given n the summation over the Dirac spec- 
trum is performed utilizing the dual-kinetic-balance (DKB) approach [|26n involving basis func- 
tions constructed from B-splines. The finite distributions of the nuclear charge and the nuclear 
magnetic moment are taken into account. 

The electric-loop part of the vacuum-polarization correction (parts (a) and (b) of Fig. ^ is 
evaluated within an effective screening potential employing the DKB method. Explicit formulas 



for the Uehling potential can be found, e.g., in Refs. [|9|, |27|]. For the evaluation of the electric- 
loop Wichmann-KroU potential we consult the approximate formulas derived in Ref. [28|- The 
magnetic-loop part (diagram in Fig. [2](c)) in the Uehling approximation was considered in Ref. 
II29I1 . We evaluate this term in the presence of an effective screening potential with account for 
the finite distribution of the nuclear charge and the nuclear magnetic moment. The value of the 
magnetic-loop Wichmann-KroU part is estimated to be rather small. It was obtained by analyzing 



its relative contribution for s states [|3 



3lh . However, in Table Hill we include its estimation to the 



total QED term for the case of 207pi,77+ 209g|78+ ^^^^ 



In the present work, the QED corrections are evaluated based on four different spherically 
symmetric binding potentials Vcs{r) which account for the interelectronic interaction between the 
valence 2pi/2 electron and the core electrons of the closed (l,s)^(2s)^ shells. The simplest choice 
of Veff is the core-Hartree (CH) potential 



POO ^ 

Kff(r) = Kuc(r) + a / dr'—p,{r') , (9) 
Jo 

where Vnuc is the nuclear potential and pc is the density of the core electrons. The screening 
potential derived from the density-functional theory reads 

Ves{r) = Kuc(r) + a / dr'—pt{r') - x^- i^rpt{r) . (10) 
Jo r> r \32iT^ J 

Here pt is the total electron density, including (ls)^(2s)^ shells and 2pi/2 electron. The parameter 
Xa is varied from to 1. The cases of = 0,2/3, and 1 correspond to the Dirac-Hartree (DH), the 
Kohn-Sham (KS), and the Dirac-Slater (DS) potentials, respectively. At large r expression (fTOl) 
should be replaced by 

v.,(o = -^^<^ (11) 

r 

to provide the proper asymptotic behavior. The self-consistent potential is generated by iterations. 

The results obtained for the QED correction for the pure nuclear potential and the screening 
potentials are presented in Table IH The first column corresponds to the values of the QED cor- 
rection for the nuclear potential V^uc- The second column gives the values of the screened QED 
correction for the core-Hartree potential Q, the other columns present the values of the screened 
QED correction for the Dirac-Hartree potential (eq. (flOl) . Xa = 0), the Kohn-Sham potential (eq. 
(fTOl) . Xa = 2/3), and the Dirac-Slater potential (eq. (flOl) . Xg = 1), respectively. The results for the 



pure nuclear potential agree with those presented in Ref . ]32l ] . 

In the calculations performed the effect of screening on QED correction is taken into account 
only in the local effective potential approximation. We estimate that for the pi/2 valence state the 
uncertainty due to this approximation amounts to about 50% for Z = 15 and decreases rapidly 
as Z increases. Evaluation of the screened QED correction within the rigorous QED approach is 
presently underway. 

The interelectronic-interaction correction of the first order in 1/Z defined by the function 
B^(aZ) can be calculated within the rigorous QED approach. Such calculations for the ground 



and first excited states of Li-like ions were performed in Refs. [|l6L 



33 



34h . The formulas derived 
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Table I: The QED correction Xj-ad for the pure nuclear potential Vnuc and for the effective screening poten- 
tials (core-Hartree, Dirac-Hartree, Kohn-Sham and Dirac-Slater potentials, respectively). 



z 




CH 


DH 


KS 


DS 


15 


0.00047 


0.00030 


0.00027 


0.00031 


0.00033 


21 


0.00041 


0.00032 


0.00029 


0.00032 


0.00033 


26 


0.00035 


0.00030 


0.00028 


0.00030 


0.00032 


37 


0.00020 


0.00020 


0.00020 


0.00020 


0.00020 


49 


-0.00005 


0.00001 


0.00002 


0.00001 


-0.00000 


57 


-0.00029 


-0.00018 


-0.00016 


-0.00019 


-0.00021 


67 


-0.00073 


-0.00056 


-0.00053 


-0.00058 


-0.00060 


75 


-0.0013 


-0.0010 


-0.0010 


-0.0011 


-0.0011 


82 


-0.0020 


-0.0017 


-0.0017 


-0.0017 


-0.0019 


83 


-0.0023 


-0.0018 


-0.0018 


-0.0019 


-0.0020 



there can easily be adopted for B-like ions by regarding the closed (Is)^ and (2s)^ shells as belong- 
ing to the vacuum state [|9|]. The results of the numerical evaluation by these formulas are presented 
in Table nil The function B'^\aZ) indicates the values obtained for the point-nucleus case using 



the method of the generalized virial relations for the Dirac equation (350. The function Bj^^^^ (ctZ) 



takes into account the nuclear charge distribution effect. Its numerical calculation is performed 
using the DKB approach. The function Bj^^\aZ) incorporates also the Bohr-Weisskopf effect. 



The higher-order term C^{Z,aZ)/Z'^ is evaluatec 



37 



using the large-scale configuration- 



3811 ■ The many-electron wave function 



interaction Dirac-Fock-Sturm (CI-DFS) method 
\l/(7 J) with the total angular momentum J and other quantum numbers 7 is expanded in terms of 
a large number of the configuration state functions (CSFs) with the same J, 



(12) 



For each relativistic atomic configuration the CSFs ^aiJ) are eigenfunctions of and can be 
obtained as linear combinations of the Slater determinants corresponding to this configuration. 
The set of the CSFs in the expansion (lT2l) was generated including all single, double and triple 
excitations. The Slater determinants are constructed from one-electron four-component Dirac 
spinors (orbitals). For the occupied shells these orbitals were obtained by the multiconfiguration 
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Table II: First-order interelectronic-interaction correction for the ground-state of B-like ions. s|°^(aZ)and 



Bj^^^\aZ) correspond to the point and extended charge nucleus, respectively. Bl^"^ ' (aZ) includes the 
BW correction. 



(BW), 



z 


Bj^\aZ) 


Br\az) 


Br\c.z) 


15 


-5.90979 


-5.9097 


-5.9097 


21 


-6.11189 


-6.1117 


-6.1116 


26 


-6.34230 


-6.3417 


-6.3413 


37 


-7.09270 


-7.090 


-7.090 


49 


-8.44727 


-8.434 


-8.431 


57 


-9.83459 


-9.801 


-9.798 


67 


-12.4859 


-12.37 


-12.35 


75 


-15.8639 


-15.56 


-15.51 


82 


-20.4837 


-19.75 


-19.51 


83 


-21.3367 


-20.49 


-20.42 



Dirac-Fock method. 

The CI-DFS method allows us to calculate the interelectronic-interaction correction to all or- 
ders in 1/Z within the Breit approximation, whereas the B^^{aZ)/Z term is obtained within the 
rigorous QED approach as described above. In order to combine these approaches, we subtract 
the value of the 1/Z term calculated within the Breit approximation from the CI-DFS result. In 
this way we obtain the C^{Z, aZ)/Z'^ contribution. 

Table Hill presents the individual contributions and the total theoretical results for the ground- 
state hyperfine splitting in B-like ions of particular interest. It can be seen that for heavy ions 
the uncertainties of the total theoretical values are completely determined by the Bohr-Weisskopf 



effect. These uncertainties can be strongly reduced employing the experimenta 
hyperfine splitting in the corresponding H-like ions. Following the works [|1(a 



values for the 



1511 we use the 



experimental values for the ground state hyperfine splitting, Ai^ixp = 1.2159(2) eV for H-like 



207t 



Pb^^+ m and AE, 



{Is) 
exp 



5.0840(8) eV for H-like ^"93^82+ extract the Bohr-Weisskopf 



corrections for the Is state employing the theoretical values for all other contributions from Ref. 
[^. Considering different models for the nuclear magnetization distribution, we have found that 
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Table III: Individual contributions to the ground-state hyperfine splitting of B-. 
error bars indicated do not include the nuclear magnetic moment uncertainties 



ike ions, in meV. The total 





455^16+ 


57pg21+ 




209gj78+ 


Effect 


= 4.7565 


= 0.090623 

^J■N 


= 0.59258 


= 4.1106 


Dirac value 


2.3126 


0.15010 


71.89 


296.35 


Finite nuclear size 


-0.0001 


-0.00001 


-2.18(1) 


-9.84(5) 


Bohr-Weisskopf 










(direct calculation) 


0.0000 


-0.00001 


-0.84(8) 


-0.97(34) 


Interelectronic interaction, 1/Z 


-0.6424 


-0.03406 


-6.82 


-28.17 


Interelectronic interaction. 










and higher orders 


0.0408(9) 


0.00182(3) 


0.24(1) 


0.98(3) 


QED 


0.0007 


0.00004 


-0.06(1) 


-0.26(3) 


Total 


1.7116(9) 


0.11788(3) 


62.23(8) 


258.09(35) 


Bohr-Weisskopf 










(from the Is experiment) 






-0.83(1) 


-1.25(4) 


Total 






62.24(2) 


257.84(5) 



(2p) /(Is) 



the ratio of the Bohr-Weisskopf corrections is rather stable, Sp{^ /e 



Pb 



0.287(2) for Pb and 

^bTV^Bi'*'* = 0.295(2) for Bi. It allows us to deduce the following values for B-like ions: Ep^^^ = 
0.0119(2) and ej^f^ = 0.00437(15). The corresponding contributions to the hyperfine splittings 
are —0.83(1) meV and —1.25(4) meV, respectively. The total theoretical values, which include 
also the modification of the Bohr-Weisskopf effect to the interelectronic-interaction correction, 
amount to 62.24(2) meV for 207pb77+ ^nd 257.84(5) meV for ^^^Bf^+. It should be stressed 
that the uncertainty of these values is not equal to the sum of the uncertainties of the individual 
contributions. This is due to the fact, that the total hyperfine splitting value found in this way is 



sufficient 



moment lIlSll . 



y stable with respect to possible variations of the nuclear charge radius and the magnetic 
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III. CONCLUSION 



In this paper we have calculated the ground-state hyperfine splitting of high-Z boronlike ions. 
The interelectronic-interaction correction is evaluated utilizing the XjZ perturbation theory and 
the large-scale configuration-interaction Dirac-Fock-Sturm method. The radiative corrections are 
calculated in the presence of an effective potential that partly accounts for the screening effect. It is 
shown that the Bohr-Weisskopf effect for heavy boronlike ions can be calculated with a high preci- 
sion utilizing the experimental value for the ground-state hyperfine splitting in the corresponding 
hydrogenlike ions. As a result, the most accurate theoretical predictions for the hyperfine splitting 
values of B-like Sc, Fe, Pb, and Bi are obtained. 
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